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Two-shell collisions in the GRB afterglow phase 
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ABSTRACT 

Strong optical and radio flares often appear in the afterglow phase of Gamma-Ray 
Bursts (GRBs). It has been proposed that colliding ultra-relativistic shells can pro- 
duce these flares. Such consecutive shells can be formed due to the variability in the 
central source of a GRB. We perform high resolution ID numerical simulations of 
late collisions between two ultra-relativistic shells in order to explore these events. 
We examine the case where a cold uniform shell collides with a self-similar Blandford 
and McKee shell in a constant density environment and consider cases with differ- 
ent Lorentz factor and energy for the uniform shell. We produce the corresponding 
on-axis light curves and emission images for the afterglow phase and examine the oc- 
currence of optical and radio flares assuming a spherical explosion and a hard-edged 
jet scenario. For our simulations we use the Adaptive Mesh Refinement version of the 
Versatile Advection Code (AMRVAC) coupled to a linear radiative transfer code to 
calculate synchrotron emission. We find steeply rising flare like behavior for small jet 
opening angles and more gradual rebrightenings for large opening angles. Synchrotron 
self-absorption is found to strongly influence the onset and shape of the radio flare. 
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1 INTRODUCTION 

The internal shock collisions version of the fireball 
model provides an adeq uate description of the origin of 
GRBs (|R ees & Mcszaros 199 4]). A coll apsing massive star 
( Wooslevlll993l;IWooslev fe Bloomll2006l ) or a binary merger 
( Naravan. Paczvnski fc Piranlll99S ) followed by a strong rel- 



ativistic explosion are considered to be the progenitors of 
these violent events. Internal collisions inside the fireball of 
relativistic shells departing from the progenitor with differ- 
ent velocities, give rise to the GRB (|Piran||2004). 

The same model attributes the afterglow emission to 
synchrotron radiation which is emitted during the decelera- 
tion of the external sho c k in t he interstellar medium (ISM ) 
i|Sari, Piran fc Naravan] 1 19981 : IZhang fc Meszaros! l2004h . 
This behaviour can last for several days or even months 
after the burst covering a wide range of the spectrum. As 
afterglow observations improved, however, certain questions 
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were r aised that cou ld not be answered with the standard 
model (|Zhand (|2007l ) , for a detailed discussio n) . Recent ob- 
servations in the o ptical Jstanek et al.l|2007|). radio as we ll 
as the X-ray band (jBurrows et al.ll2005l ; lNousek et alj|2006h . 
show a strong variability in the afterglow phase for a large 
proportion of the bursts, which can not be reproduced by 
the standard external shock model. 

It has been proposed that a bump in the afterglow light 
curve may result when the forward shock propagating in 
the ISM encounters a density jump, caused by an inhomo- 
geneity of the surrounding medium generated by interstellar 
turbulence or by anisotropy in a pre c ursor wind from th e 
GRB progenitor l|Wang fc Loebl l2000l ; lEazzati et al]|2002T l. 
Bowever, numerical simulations of a spherical explosion ex- 
hibit a rather canonical behaviour and even for a sharp 
and large increase in the external density this model does 
not produce sharp features in the light curve and cannot 
accou nt for significant temporal variability in GRB after- 
glows (|Nakar fc Granotj[2 006 ; van Eerten et al1l2009l ). It has 
been suggested that a late activity o f the central engine 
could explain the observed variability llFalcone et al.ll2006l : 
IZhang et al]|2006l ; iRomano et alll2006l : lloka et al.ll2005l ). 
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In the late activity scenario, the central engine produces 
consecutive explosions after the initial burst which collide 
when a slow shell is followed by a faster one. This late ac- 
tivity of the source could be explained from a two-stage 
collap se in the central object. As proposed by iKing et al.l 
( 2005) a collapsing core which has enough angular momen- 
tum to fragment will leave behind a second compact "star" 
in the form of a self-gravitating neutron lump. The fallback 
of this "star" at later times on the initial compact object can 
restart the central engine. Other theories suggest that the 
viscous hyperaccreting accretion disk around a black hole 
which fragments at large radii becomes dynamically unsta- 
ble on different timescales and thus c ollapses at different 
times (|Perna. Armitage fc Zh ang 2006). It is proposed that 
the region at the vicinity of the accretor can play an impor- 
tant role in determining the acc retion rate and therefo re the 
energy output of the explosion |Proga fc Zhangj [2006'). 

According to late activity models the second blast wave 
continuously supplies the system with energy while colliding 
with the initially ejected material, producing in that way the 
rebrightening observed in the afterglow. 

The role of magnetic fields in GRBs is still arguable. 
In the fireball model the presence of a magnetic field is 
not dynamically important for the evolution of the flow 
but plays an important role for the emission during the in- 
teraction of the flow with the external medium. Although 
the early afterglow emission strongly depends on the mag- 
netization of the flow, in the late stages of the afterglow 
where the shells experience strong deceleration, the evo- 
lution of strongly magnetized shells resembles that of hy- 
drodynamic shells and can be described by the se lf-similar 
Blandfor d-Mckee (BM) dBlandford fc McKeeil 19761 ) approx- 
imation (|Mimica et al.ll2009h . At this stage of the afterglow 
the emission no longer contains information about the initial 
magnetization of the flow. 

In section 2 we describe the high resolution numeri- 
cal simulations we performed of late collisions of two ultra- 
relativistic shells during the afterglow phase. We claim that 
differences in the flow must have an impact on the result- 
ing light curves and perform four simulations with varying 
Lorentz factor and energy content of the second shell in order 
to investigate the effect of these parameters. The adaptive 
mesh refinement (AMR) technique enables us to use high 
resolution in long term, ID relativistic hydro simulations, in 
order to capture the forward and reverse shock formation 
on the second shell and study in detail the stages before, 
during and after the merger of the two shells. 

The effects of the collision between the two shells in the 
light curves is described in section 3. We study both spheri- 
cal explosions as well as a hard-edged jet scenario where no 
lateral spreading has occured (numerical simulations in two 
dimensions have shown that only very modest lateral spread- 
ing occurs while the jet is relativistic dZhang fc MacFadvenl 
120091 ; iGranot et aLll200ll ; iMeliani fc Keppenj|2010h . Qpti- 
cal and radio on-axis light curves are calculated for different 
opening angles and the strength of the occuring flare or re- 
brightening is found to depend on this opening angle. We 
also note a clear difference in the shape between optical and 
radio light curves as well as a difference in the time of the 
appearance of the flare between the two frequencies. We ex- 
plain this chromatic behaviour in terms of the synchrotron 
self-absorption mechanism and the different main contribut- 



ing regions of the jet to the emission. We construct emission 
images for the different stages of the merger and connect the 
dynamical characteristics of the flow at each stage of the col- 
lision to the features in the light curves. We will discuss and 
summarize our results in section 4. 

For the dynamical simulations we are using the Adap- 
tive Mesh Refinement version of the Ve r satile Advec- 
tion Code (AMRVAC) l|Keppens et al.ll2003l ;l Meliani et al.l 
l2007h . and for the light curves and emission images calcula- 
tions the radiation code of van Eerten & Wijers (2009). 



2 MODELING OF THE MULTI-SHELL 
DYNAMICS 

When the initially ultra-relativistic shell ejected from the 
central source starts to decelerate in the interstellar medium, 
a forward shock is created separating the shocked ISM from 
the ambient ISM. As mass is swept up, the kinetic energy 
of the shell is transformed into kinetic and thermal energy 
of the shocked matter. At the same time a reverse shock is 
formed which crosses the shell leading to conversion of the 
shell's kinetic energy into thermal. The resulting shocked 
ISM matter ultimately follows the self-similar BM analyti- 
cal solution (|Meliani et al.|[2007l 'l . The afterglow is nowadays 
widely recognized as synchrotron radiati on emitted during 
this p hase of the propagation of the shell (jMeszaros fc Reea 
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7). 



In our model we consider that the central engine re- 
mains active even after the initial ejection of the first shell 
resulting in a delayed second explosion. The produced blast 
wave will now travel with a steady velocity into an empty 
medium, since most of the matter has been swept up, until 
it reaches the termination shock of the first shell. In this 
paper we reproduce the collision process of these two shells 
and claim that for small opening angles the heating of the 
matter that happens during this phase is responsible for the 
appearance of the flares observed in the light curves. 



2.1 Special relativistic hydrodynamic equations 

We perform the dynamical simulations using the ID special 
relativistic hydrodynamic equations in spherical coordinates 
and the code AMRVAC. The equations describing the mo- 
tion of a relativistic fluid are given by the five conservation 
laws 



(pu"),„ = 0, (TH; 







(1) 



where p,,u = 0,1,2,3 are the indices running over the 4- 
dimensional spacetime, p is the proper rest mass density 
of the fluid, u M is the four-velocity and is the stress- 
energy tensor given by = phu^u v + pg 1 *". Here with 
p we denote the fluid rest frame pressure, while g^ v is the 
Minkowski metric tensor, as we will consider a flat space- 
time at distances far from the central engine. The specific 
enthalpy h of the fluid is given byh=l+e+p/p where 
e is the specific internal energy. Rewriting the conservation 
equations in vector form we have in a familiar 3+1 split the 
conservation laws 



&U_ dF l {U) 
~dt + dx* 



= 0, with i = 1,2,3. 



(2) 
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The vector U is defined by the conserved variables as 
U = [D = pj, S = ph-y 2 v, r = ph-y 2 — p — D] T , 
and the fluxes are given by 



F = [pyv, p/17 2 vv + pi, ph-y 2 v - 



■pjv] 



(3) 



(4) 



where v is the three-velocity and I is the 3x3 identity ma- 
trix. The system of equations is closed by using the equation 
of state 



p = (T P - \)pe. 



(•>) 



In our simulations we choose the polytropic index to be T p = 
4/3. This is a good approximation since most of the shocks 
in the cases described below are mainly relativistic or near- 
relativistic. 



2.2 Initial setting 

For the dynamics of the first shell we consider that the re- 
verse shock has already crossed the shell which now deceler- 
ates in the external medium. For the purpose of this simula- 
tion we will use the Blandford & McKee (BM) approxima- 
tion to describe this phase. This is a self-similar solution of 
a relativistic blast wave expanding in a uniform or radially 
varying medium. We consider the case where the explosion 
is assumed to be spherically symmetric and adiabatic. 

In the BM model the density of the circumburst medium 
scales as a power law with distance pi(r) oc r~ k . For 
all the simulations in the present paper we will consider 
that the density of the circumburst medium is constant 
(k — 0), with particle number density m = 1cm -3 and 
cold, with the pressure given by pi — W~ 5 nirn p c 2 chosen 
such that it does not dynamically affect the system. We 
set the Lorentz factor of the BM shock at To = 23 at the 
start of the simulation placed in distance Ro — 2.04 x 10 17 
cm. Considering the decelerating radius of the BM shock 
Rdec = (3-Eo/(47rnom p c 2 r 2 )) 1 ' /3 , after which the Lorentz 
factor of the shock starts decreasing with distance as a power 
law, the initial distance Ro of the shock corresponds to a dis- 
tance 3.7 times greater than the deceleration radius Rdec of 
a jet with initial Lorentz factor 100 and 7.8 times greater 
than the deceleration radius of a jet with initial Lorentz fac- 
tor 300. The energy content of the shell is E — 10 52 erg. 
According to the BM model in the ultra-relativistic case the 
jump conditions at the BM shock are given by 

2 o 
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where index 2 denotes the shocked medium and V is the 
Lorentz factor of the shock. According to the BM model in 
the ultra-relativistic case, the radius of the shock at time t 
is given up to order 0(F -2 ) by 

1 

8~r~ 2 



R (t) = C M 1 - ^2 



(9) 



From the jump conditions and by choosing the similarity 
variable to be x = [l + 8F 2 ] (1 - r/t), we obtain the prop- 
erties of the shocked medium 



Table 1. Properties of the second shell for each case. T = 23 and 
E = 10 52 ergs are the Lorentz factor and energy of the BM shell. 
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The total energy is then given by E = 87rpiFgc 5 ^/17. If 
the initial Lorentz factor of the simulation is fixed at To, 
the duration to of the shock so far then follows from this 
equation. The initial pressure and density jumps between 
the BM shell and the ISM at the position of the shock are 
P2/P1 — 10 7 and P2/P1 = 10 2 . 

The second shell is uniform, cold and ultra-relativistic 
and placed at distance AR = 10 14 cm behind the BM shell. 
It is therefore assumed that the shell has moved freely up 
to this point. Considering a duration of the second ejection 
event of At — 1000 s, the initial thickness of the shell will 
be S = cAt — 3 x 10 13 cm. The energy of the second shell is 
given by 



E sh = ^T 2 3h R 2 n 5p sh c , 



(13) 



where Ri n denotes the initial distance of the shell and p s h 
and Y a h the initial density and Lorentz factor. The initial 
pressure is chosen as p s h = 5 x 10~ 2 p s h- In our initial con- 
ditions, we vary from case to case the given parameters 
and E s h which then serve to specify the shell density p s h. 
The energy provided here refers to the isotropic-equivalent 
energy. For all cases described in this paper we consider 
emission along the rotation axis of the system and that the 
two ejecta have the same opening angle. We also neglect 
the effects of lateral spreading in both dynamical and radia - 
tion calculations. As shown in IZhang fc MacFadvenl (|2009h . 
sideways expansion can be a very slow process and definitely 
negligible for times under consideration in this paper. 

We perform four simulations with varying Lorentz fac- 
tor and energy for the second shell. In this simulation we are 
using a domain of size [0.01, 10] x 10 18 cm and 240 cells at the 
coarsest level of refinement. The physical properties of the 
afterglow shock collision model require a large domain and 
a very thin second shell which demands very high resolution 
in order to be resolved. The maximum level of refinement is 
22, leading to an effective resolution of 5.03316 x 10 8 cells. 
We force the front and the back of the uniform shell to be 
refined in order to avoid numerical diffusion which would 
cause an artificial spreading of the shell. A summary of the 
simulation parameters can be found in Table 1. The top left 
figures in Figs. Q] and [2] show the initial conditions of case 1 
and 4 respectively. 



4 A. Vlasis, H. van Eerten, Z. Meliani, R. Keppens 



2.3 Interaction between two shells. 

2.3.1 Jump conditions. 

The evolution of the two-shell system is shown in figures 
[T] and [2] for case f and case 4 respectively. In case 1 the 
two shells initially have the same energy and Lorentz factor. 
As the BM shell decelerates in the interstellar medium the 
second shell catches up. While the matter from the BM shell 
is swept up by the second shell, a forward shock is created 
separating the shocked matter from the unshocked. At the 
same time a reverse shock crosses the second shell and a 
contact discontinuity appears in between both shocks. At 
this stage the front of the second shell has split into four 
regions. In region 1 there is the BM matter. In region 2 there 
is the BM matter which has been heated by the forward 
shock. In region 3 there is the matter of the uniform shell 
that has been heated from the propagation of the reverse 
shock and in region 4 there is the unshocked matter of the 
second shell. This is also shown for all 4 cases in Fig. [3] where 
a close-up of the second shell at early times is given. 

The jump cond it ions at the second shell as given by 
iBlandford fc McKed l| 19761 ) for an arbitrary strong shock 
{p2/n 2 > Pi/ni) are 

e 3 =7 a — — =4t 2 +3, (14) 
ni ni 

-^=73 -1, -= 47 3 + 3, (15) 

where 7 2 is the Lorentz factor of the fluid in region 2 relative 
to region 1 and 7 3 is the Lorentz factor of the fluid in region 
3 relative to region 4. The primitive variables iii , pi as well as 
the internal energy density d and enthalpy uii , are measured 
in the fluid frame, while the Lorentz factors of the several 
regions, 7 y i , are measured in respect to the ISM which is 
considered to be at rest. In our case the forward shock of the 
second shell is moving in a hot medium which has already 
been heated by the BM shock. Therefore wi = ei +pi — 4pi, 
assuming the ultra-relativistic equation of state ei ~ 3pi. 

2.3.2 Dynamics of the different cases. 

After the initial explosion and acceleration of the uniform 
shell, the swept up BM matter starts playing an effective role 
in the deceleration of this shell. In Fig.[3]we show snapshots 
of all 4 simulations taken at emission time t e = 6.98 x 10 6 s 
after the explosion when the forward shock, the contact dis- 
continuity and the reverse shock are fully developed. At this 
time the BM matter is continuously heated by the forward 
shock resulting in an increase in the density of the shocked 
matter, n^/ni ~ 7.10 (case 1). The relative Lorentz factor 

is given by 7 2 = 7271(1 - V( X ^7 2 ~ 2 )( 1 ~ 7i~ 2 )) ~ L46 - A 
small inconsistency is observed between the simulation re- 
sults and the jump conditions which derives from the fact 
that in all our simulations the random kinetic energy per 
particle at the two sides of the shock is P2/712 > pi/ni rather 
than P2/712 2> pi/ni. Similarly the propagation of the for- 
ward shock results in an increase of the internal energy of the 
shocked matter which now is e 2 ~ 0.61 ~ ^j 2 w\n-ijn\. At the 
same time the reverse shock crosses the shell while heating 
the uniform shell matter and transforming the kinetic energy 
of the shell into thermal. The contact discontinuity separat- 
ing the two regions appears as a density jump between the 



two shocked fluids, 713/712 — 10 2 while the Lorentz factor 
and pressure remain continuous, 72 = 73 and P2 = P3- The 
efficiency of this energy transformation mechanism is highly 
dependent on the reverse shock itself. According to Sari & 
Piran (1995) the reverse shock depends only on two parame- 
ters, the Lorentz factor of the unshocked shell material rela- 
tive to the BM unshocked matter, 7 4 , and the density jump, 
714/721. For 74 3> 714/ni the reverse shock will be relativistic, 
7 3 S> 1. For case 1, at emission time t e = 6.98 x 10 6 s. (Fig. 

the reverse shock remains Newtonian (73 = 1.00484) 
while propagating into the shell and thus insufficient to heat 
the shell effectively. Following this behaviour, by the time 
the reverse shock reaches the back of the shell, the shell is 
still dominated by its kinetic energy although significantly 
decelerated (see also Fig. [TJ. 

We follow the same analysis for case 2 in which we dou- 
ble the Lorentz factor of the second shell compared to case 
1. Since we choose to maintain constant the energy and the 
thickness of the second shell, the change in the Lorentz factor 
affects the density of the second shell, which is now smaller 
by a factor of 4 compared to case 1 (eq. U3[) . Comparing the 
properties of the flow for the two cases at similar times after 
the explosion we notice that the forward shock is stronger in 
case 2. Specifically at emission time t e = 6.98 x 10 6 s after 
the explosion the forward shock has a Lorentz factor rela- 
tive to the BM matter 7 2 = 1.99 and the number density 
satisfies 712/711 ~ 7.67. That shows that the forward shock is 
more efficient in the second case as it compresses the matter 
of the BM shell to a higher degree and to higher Lorentz 
factor than in case 1. At the same time the reverse shock 
is stronger than in case 1, 7 3 = 1.16 and although it is not 
ultrarelativistic it is more efficient in converting the kinetic 
energy of the second shell into thermal. That appears clearly 
as a difference in the pressure between the shocked and un- 
shocked shell matter, pzIpa = 12.27 for case 2 compared to 
Ps/Pi = 1-78 for case 1. 

In case 3 we double the energy of the second shell. The 
relevant Lorentz factor between the shocked and unshocked 
BM matter is 7 2 ~ 1.75 and the density jump satisfies 
712/711 ~ 6.69. Compared to case 1 the reverse shock for 
this case remains very weak, 73 ~ 1.00013, while propa- 
gating in the uniform shell since it has to traverse a denser 
medium than before, leading to a compression ratio between 
the shocked and the unshocked shell matter PzIpa = 1.09. In 
case 4 we double both the energy and Lorentz factor of the 
second shell. The initially very fast shell leads to a strong 
forward shock, 7 2 ~ 2.26. The reverse shock is now stronger 
compared to case 3, 73 — 1.15, compressing the shocked 
shell matter to higher pressure, pi, I Pa = 7.98. 

In figure [3] we plot the thermal to mass energy ratio 
together with Lorentz factor, density and pressure for fixed 
early time for all cases. We observe that as matter from the 
BM shell is swept up from the forward shock of the second 
shell, the thermal energy of the fluid increases compared 
to the mass energy (which appears as a small bump in the 
thermal to mass energy curve at the position of the forward 
shock). While the shell gets traversed by the reverse shock 
part of the kinetic energy of the shell is transformed into 
thermal as a result of the propagation of the reverse shock. 
We notice that the relativistic reverse shock in case 2 is a lot 
more efficient in raising the thermal energy component than 
in case 1. In case 3 and 4 a similar behaviour is observed. By 
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Figure 1. Snapshots of the dynamics for case 1, taken at local emission times t e = 6.81 X 10 6 s, t e = 8.48 X 10 6 s, t e = 1.15 X 10 7 s and 
t e = 2.35 X 10 7 s after the initial explosion (top left - bottom right). Lorentz factor, density and pressure are indicated as the dotted, 
solid, and dashed line. The top-left figure indicates the initial setting of our simulation. 



increasing the energy in case 3 and with the Lorentz factor 
being the same as in case 1, the density becomes higher in 
the shell (eq. 13). As a result the reverse shock is very weak 
and highly inefficient in thermalizing the cold shell. In case 
4 both the energy and the Lorentz factor of the second shell 
have twice the value compared to case 1. In this case the 
reverse shock is propagating in a less dense medium com- 
pared to case 3 and is now more efficient in transforming 
the kinetic energy of the second shell into thermal (Fig. [3]). 
The Lorentz factor of the reverse shock is slightly overesti- 
mated as a result of our use of a fixed adiabatic index of 4/3 
throughout the domain. Nevertheless, as we see from Fig. O 
the fluid is highly relativistic from the contact discontinuity 
of the second shell onward and therefore the regions that 
mostly contribute to the observed flux (see next section) are 
appropriately described by a fixed adiabatic index. 

Later, at emission time t e = 1.513 X 10 7 s (see Figs. [1] 
and [2]), the forward shock will catch up and overcome the 
BM shell as the latter one decelerates in the ISM. As the 
forward shocks merge and the reverse shock has crossed the 
back of the shell, a dense but slow and underpressured region 
is left behind unable to follow the forward shock as it has lost 
almost all of its kinetic energy reaching a near-equilibrium 
state with the surrounding matter. 



3 RADIATION CALCULATIONS 

In this section we describe the numerical calculations we per- 
formed in order to construct the light curves at the afterglow 
phase. The following calculati ons were carried out wit h the 
radiation code introdu ced in Ivan Eerten et all (|2009l ) and 
Ivan Eerten et al. (2010). The main process contributing to 
the afterglow is synchrotron radiation. In both external and 
internal shock collision models, a magnetic field is required 
in order to fit the observational data. This magnetic field is 
most likely generated by instabilities forming in the shock 
such as the relativistic two-stream instability or th e Weibel 
instability jMedvedev fc Loeb|[l999l ; IWeibelll959h . In this 
case a fraction of the total thermal energy behind the shock 
goes to particle acceleration and another one to the gen- 
eration of the magnetic field. Assuming that the fractions 
of the thermal energy density e t h, contributing to the mag- 
netic energy density, es, and electron energy density, e_e, 
have a fixed value, we can calculate the afterglow emission. 
In our calculations we neglect the effect of the magnetic 
fields and radiative losses on the dynamics and do not take 
into account the effects of Compton scattering and electron 
cooling (that play a negligible role at times and frequencies 
under consideration). We do however include the effect of 
synchrotron self-absorption (ssa) due to the re-absorption 
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Figure 2. Snapshots of the dynamics for case 4, taken at local emission times t e = 6.81 X 10 6 s, t e = 8.48 X 10 6 s, t e = 1.15 X 10 7 s and 
t e = 2.35 X 10 7 s after the initial explosion (top left - bottom right). Lorentz factor, density and pressure are indicated as the dotted, 
solid, and dashed line. The top-left figure indicates the initial setting of our simulation. 



of the radiation from the synchrotron electrons by solving 
simultaneously the linear radiative transfer equations for a 
large number of rays through the evolving fluid. An analyt- 
ical formula for obtaining the self-absorption frequency v sa 
can be found in lGranot. Piran fc Saril (|l999bT ) . 

The radiation code is specifically written to include the 
snapshots produced by the dynamical simulations performed 
with AMRVAC. In all our calculations the values of es and 
€e are fixed to 0.01 and 0.1 respectively and we assume 
a power law distribution for the accelerated electrons with 
p = 2.5. We also fix the fraction of the electrons accelerated 
to this power law distribution £jv equal to 0.1. Assuming 
isotropic radiation in the comoving frame, the observer's 
flux calculated by the radiation code is given by 



F„ 



l + z 



d 2 P^ 



dudQ. 



(1 - Pn)cdAdt e 



(16) 



in the optically thin limit, where for the purpose of our sim- 
ulations the redshift z is chosen to be zero. Here Dl denotes 
the observer luminosity distance, f3 the fluid velocity in units 
of c, d 2 Py I dudQ, is the received power per unit volume, fre- 
quency and solid angle and fi — cos 9, with 6 being the angle 
between the fluid velocity and the line of sight. The surface 
element dA corresponds to an equidistant surface A, which is 
a surface intersecting the fluid grid from which the radiation 



arrives at the observer at time t„b s . Each surface corresponds 
to a specific emission time t e . For a photon emitted from a 
location (r,#) at emission time t e , the observer time is given 
by 



tobs 



t e 



Vfl 

c 



(17) 



When electron cooling does not play a role, the shape 
of the observed spectrum follows directly from the dimen- 
sionless function Q[y/v m ) which has the limiting behaviour 
Q <x (v/v m ) 1/3 for small [v/v m ) and Q tx (y /v m ) i - 1 ^ v)/2 for 
large (y/v m )- The received power depends on this shape and 
on the local fluid quantities via 



d 2 P v ^ N nB' 
dudQ. K 7 3 (1 -P/i) 3 



(18) 



where n is the lab frame number density of the electrons and 
B' is the comoving magnetic field strength. The synchrotron 
peak frequency u m corresponds to the Lorentz factor of the 
lower cut-off of the accelerated electron's power-law distri- 
bution 7 m , assuming that the Lorentz factor for the upper 
cut-off goes to infinity. Then the lower cut-off for the elec- 
trons will relate to the comoving number density ri and 
thermal energy density e' th via 
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Figure 3. Ratio of thermal to mass energy (E t f l /(pc 2 )) and normalized Lorentz factor, density and pressure for case 1 to case 4 (top 
left to bottom right) for the second shell at emission time t e = 6.98 X 10 6 s. Distance is normalized to 10 17 cm. For cases 1 and 2 the 
transformation from thermal to kinetic energy is clear at the forward shock as well as the kinetic to thermal energy transformation at the 
position of the reverse shock. The ratio depends strongly on the properties of the second shell. The four regions of interaction between 
the two shells are indicated on the figure for case 1. Region (1) consists of the unshocked BM matter, region (2) of the shocked BM 
matter, region (3) of the shocked uniform shell matter and region (4) of the unshocked uniform shell matter. We focus on the front of 
the second shell, the dynamical evolution of which plays an important role on the light curves. 



7m OC j. (19) 

\P - 1/ tivn 

The particle distribution and its lower cut-off are set 
at the shock front. However, subsequent evolution of 7 m 
is dictated by adiabatic expansion of the fluid rather than 
synchrotron radiative losses. Therefore, in a relativistic fluid, 
eq. (|19|l also holds further from the shock front with the same 
value of tE- 

The temporal behaviour of several GRBs contradicts 
the spherical explosion senario. A rapidly decaying afterglow 
emission suggests that the flow must be collimated rather 
than spherical. This is vital for the GRB mechanism since 
a spherical expansion would require a total energy budget 
~ 10 54 ergs which is hard to produce from a stellar mass 
progenitor, while a jet shaped explosion can have the same 
result with less energy, ~ 10 51 — 10 52 ergs. This model sug- 
gests that when the jet decelerates to Lorentz factors such 
that 7 ~ O^ 1 is satisfied, with 8h being the half jet opening 
angle, the flux that the observer receives will start decreas- 
ing resulting to a break in the afterglow light curves. In 



contrast to AGN jets that can be directly observed, GRB 
jets are only implicitly assumed from this break. 

In our simulations we construct optical and radio light 
curves and emission images at various times for the 4 cases 
described above. We try different opening angles of the jet 
and associate the flare characteristics with the jet opening 
angle. 

The main characteristic of our approach is the sepa- 
ration of the dynamics from the radiation simulation. The 
outputs from the dynamical simulations are used as an input 
to the radiation code in order to calculate synchrotron emis- 
sion. When absorption plays a role, the code solves a series 
of linear radiative transfer equations (rather than eq. 16 di- 
rectly) along the light rays starting from the back of the jet 
and passing through the jet towards the observer. For given 
emission time t e a surface on the radiative volume of the jet 
exists from which emission arrives at the observer at time 
t bs- A summation over the light rays along this equidistant 
surface, and an integration over all the equidistant surfaces 
is applied in order to calculate the contribution to the emis- 
sion of every part of the jet. The observed flux is obtained by 
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summing over the rays emerging from the jet. When the rays 
are not summed over, a spatially resolved emission image 
of the two shell system is obtained. An adaptive procedure 
similar to the one used in AMRVAC is employed when more 
rays need to be calculated in order to adequately capture 
the emission from the underlying fluid profile. More details 
on the radiation c ode algorith mic strategy can be found in 
Ivan Eerten et all |2009l . l2010l ). 



Table 2. Observer time estimation of the jet break (estimation 
made in days after the explosion) for the four simulated cases 
assuming two hard-edged jet scenarios. 

Jet break estimation 





case 1 


case 2 


case 3 


case 4 




= 2 


0.0028 


0.0028 


0.0035 


0.0035 


2d h 


= 5 


0.032 


0.032 


0.041 


0.041 



3.1 Early afterglow and jet break estimation. 

The dynamical simulations cover a timescale starting from 
0.072 days and ending 10 days after the initial explosion 
in the observer's time frame. In order to include the initial 
stages of the afterglow in our simulation we assume that 
prior to the simulation the outer shock has evolved according 
to BM and the second has moved with a constant velocity 
while retaining its initial shape. In this way the deceleration 
of the initial explosion is taken into account resulting in the 
appearance of the jet break in the light curve. 

Due to relativistic beaming the emission that reaches 
the observer is limited to emission angles 9 < 7 -1 . During 
the afterglow phase however the flow is significantly deceler- 
ated and thus larger emission angles can be observed. When 
the Lorentz factor becomes small enough so that the con- 
dition 6h ~ 7" 1 is satisfied, where 6h is the half opening 
angle of the jet, the observation cone becomes big enough 
for the edges of the jet to become visible to the observer. 
When there is no significant spreading of the matter at the 
edges only part of the visible region is occupied by the jet 
and thus from this point after the flux will start decaying 
faster. This transition to the faster decaying part of the light 
curve is often referred to as the je t break. 

An analytical formula given in lvan Eerten et ail |2010bl ) 
gives an estimation of the observed jet break time, t bs,br 
depending on the total energy of the explosion E, the half 
opening angle 9h of the jet, the circumburst density no and 
profile k. The contributing area to the emission in this case 
is not only the shock front itself but also the area behind 
the shock. For that purpose the radial profile from the BM 
model has been used, and the formula reads 



, _ fl 2+2/(3-fc) 

l>obs,br — 



1 

- + 



At 



2(1 + 2(4 -fc)) 



+ 



V(3-fe) 



2 2(1 + 2(4 -k)) 2(4 - Jfe)2(l + 2(4 - k)) 



(20) 



where k corresponds to the parameter defining the power- 
law of the circumburst medium density, which for our case is 
always set to 0, and A\ a parameter of the fluid depending 
on the energy E of the explosion and the number density at 
the position of the shock, Ai = £(17-4fc)/(87rm p no7?§c 5_fc ) 
with no being the proton number density at distance Rq. 

The jet break observer time estimation t b s ,br, is pre- 
sented in Table [2] for all the cases we simulate and for two 
different opening angles of the jet. As expected, when the 
total energy of the explosion remains the same, the jet break 
time remains the same even if the Lorentz factor of the sec- 
ond shell is different. For larger opening angles of the jet 
the jet break is observed at later observer times as it takes 
longer for the shell to decelerate to Lorentz factors such that 



Oh ~ condition is satisfied. We discuss the jet break 
characteristics further in section 3.4. 

In our model the interaction between the two shells 
starts happening after the jet break has occured, leading 
to a rebrightening of the afterglow and a sudden increase in 
the flux. The shape of the light curve and the characteristics 
of the emission are subjected to the dynamics of the flow and 
present qualitative and quantitative differences from case to 
case. 



3.2 Optical and radio light curves 

For the four cases described in section 2.3 we produce optical 
and radio light curves. We examine each case assuming dif- 
ferent opening angles of the jet. In Figure [4] and [5] we present 
the light curves produced by the radiation code for spherical 
and collimated expansion of the two colliding shells. 

A rebrightening of the optical light curve is observed for 
tobs = 0.2 days for the spherical explosion in case 1 (Fig. 3}, 
which is attributed to the interaction of the two shells when 
the second one reaches the back of the BM shell. We find 
that the smaller the opening angle of the jet, the steeper the 
rise of the flare. This can be understood from the fact that 
for a jet with small opening angle there is, at a given ob- 
server time, less emission still arriving from earlier emission 
times and higher emission angles. As we show in section 3.5 
(Figures [9l I10p early emission time contribution comes from 
high emission angles on the jet, thus the features of the light 
curve are less smeared out for a collimated outflow. 

This can also explain the difference we observe in the 
decreasing rate of the flux at late observer times. Since there 
is no contribution from early emission times, the flux for a 
jetted outflow will decline faster compared to a spherical 
explosion. We also observe that the flare is sharper in case 2 
compared to case 1 and in case 4 compared to all four cases 
(Fig. [4] and [5]). This confirms our hypothesis that the flare 
is strongly dependent on the dynamical properties of the 
collision. For case 2, where the forward and reverse shock of 
the second shell are significantly stronger compared to case 
1 (Fig. [3} , the change in the flux during the collision of the 
two shells appears substantially stronger. In cases 3 and 4 
we observe the same behaviour in the light curves. In case 
4, the higher Lorentz factor together with the higher energy 
imposed in the second shell, lead to a flare clearly stronger 
compared to all the other cases for a 2 degrees opening angle 
of the jet. In that case the timescale of the variation for the 
flare is AT/T = 1.08 for the optical and AT/T = 1.90 for 
the radio, where AT is calculated as the full width at half 
maximum of the flare and T is the observer time of the peak. 
The relative flux increase with respect to the underlying 
afterglow is AF/F = 3.95 for the optical and AF/F = 
1.04 for the radio which is significantly reduced due to ssa 
mechanism. Although extensive comparison to observational 
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Figure 4. Optical (top) and radio (bottom) light curves for case 1 (left) and case 2 (right) for different values of the jet opening angle 
29^. In all cases an increase of the flux due to energy injection from the second shell is observed. For small opening angles a flare appears 
in all four simulated cases which differs in shape according to the frequency. The peak flux is greater the higher the Lorentz factor and 
energy of the second shell are. The jet break estimation is denoted with a vertical line for both 2 and 5 degree jets. For case I we overplot 
the light curves produced for a single BM shell for a spherical explosion and a 2 degrees hard-edged jet. The analytical estimation of the 
slope is shown above the optical and radio light curves. 



data is beyond the scope of this work, we note that this 
case resembles the afterglow of GRB 060206 which shows 
an increase in brightness by ~ 1 mag 1 hr af ter the burst, 
followed by a typical broken power-law decay (|Stanek et al.l 
120071 : IWozniakfl2006bh . 

Before and after the flares the resulting light curves fol- 
low the analytically predicted slopes for a single forward 
shock as shown in Fig [4] The light curves at a given fre- 
quency depend on the temporal evolution of the charac- 
teristic frequencies of the system, that is in our case the 
synchrotron peak frequency v m and the self-a b sorpt ion fre- 
quency v sa . As described in iGranot fc Sarll (|2002T l these 
slopes read t 1 ^ 2 for times before the break frequency u m , 
crosses the observed frequency v and t 3 ( 1-p )/ 4 for the time 
after. In addition to this behaviour the slope of the light 
curve steepens by t -3//4 after the jet break, in both opti- 
cal and radio frequency, for small opening angles of the jet. 
This steepening appears clearer in the radio (t t s ,br = 0.0054 
days for 2 deg. jet and t b s ,br = 0.068 days for 5 deg. jet) 
where the emission comes from late emission times and lower 
emission angles close to the jet axis. In the optical, where 
earlier emission times also contribute to the observed flux 



(as discussed in section 3.5 and shown in Fig. [9l and [TO)) . 
this steepening is delayed and appears at t b s = 0.075 days 
almost at the same time as the break frequency v m crosses 
the optical band. 

Throughout our calculations and for all cases the critical 
frequencies satisfy v c > v a b s , v c > v m and v c S> v sa at the 
BM shock and the forward shock of the second shell, where 
v c is the cooling frequency (all critical freque ncies are calcu- 
lated using the formulas found in Table 2 of lGranot fc Sari 
( 2002) for a BM solution. Between t obs of 0.001 and 10 days, 
v c decreases from 10 20 Hz to 10 17 Hz, while v m decreases 
from 10 16 Hz to 10 10 Hz for the BM shock of case 1). Hence, 
for the frequencies discussed in the present work we are al- 
ways in the slow cooling regime. In addition to the fact that 
the shock regions forming at the second shell are a lot thin- 
ner than a typical BM profile, this allows us to neglect, with 
small error, the effect of synchrotron cooling on the electron 
distribution for the times under consideration and conse- 
quently any changes on the self-absorption coefficient. 

Furthermore our assumption that inverse Compton (IC) 
doesn't influence the observed spectrum is confirmed as fol- 
lows. There are two ways in which IC scattering can change 
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Figure 5. Optical (top) and radio (bottom) light curves for case 3 (left) and case 4 (right) for different values of the jet opening angle 
29 n . In all cases an increase of the flux due to energy injection from the second shell is observed. For small opening angles a flare appears 
in all four simulated cases which differs in shape according to the frequency. The peak flux is greater the higher the Lorentz factor and 
energy of the second shell are. The jet break estimation is denoted with a vertical line for both 2 and 5 degree jets. 



the overall synchrotron spectral component. First by pro- 
ducing an additional emission component at high frequen- 
cies and second by dominating the electron cooling and 
thus reducing the av ailable energ y for s ynchrotron radia- 
tion. As discussed in I Sari fc Esinl (|200ll ). the first is esti- 
mated by considering the ratio of specific fluxes measured 
at the peak of the respective flux components, f m ax for the 
synchrotron and f^ ax for the IC respectively. Then this ra- 
tio is fmax/ fmax ~ l/SarnR, where ctt is the Thompson 
cross-section and n the electron density at distance R. This 
ratio remains well below unity for all the shock surfaces 
in all four simulations which are taken into consideration. 
The second factor can be estimated directly from the ratio 
rj^E /es ■ As shown by Sari & Esin, the IC cooling rate will be 
unimportant compared to synchrotron if rjeE /es <C 1, where 
V = (7c/7m) 2_p = {vc/vm)~^ p ~ 2 ^ 2 ■ In our case this ratio is 
close to unity for the early afterglow but decreases rapidly 
and remains below unity during our simulation covering the 
pre and post flaring activity period. 

3.3 The shape of the optical flare. 

The flaring activity as observed in the optical follows three 
stages. At the first stage the Lorentz factor and the thermal 
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Figure 6. Optical light curve and single emission images for sev- 
eral snapshots of the two shell system for case 4. The emission 
time of the snapshots range from t e = 7.64 X 10 6 s for snapshot 
50 to t e = 1.26 X 10 7 s for snapshot 350. 
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energy of the BM matter in front of the second shell increase 
due to the propagation of the forward shock leading to the 
sudden rise of the flux. The flare observed in Fig. [4] and is 
attributed to that increase. This behaviour continues until 
tots = 0.23 days when the reverse shock crosses the back 
of the second shell. At that time the forward shock starts 
decelerating and that can be seen as a reduction of the flux, 
after the initial peak, in the optical light curves. At time 
tobs = 0.25 days the shocked BM matter has separated from 
the shell and now propagates in the BM shell while con- 
stantly heating the swept up matter. The flux for that stage 
of the motion remains almost constant in the optical light 
curve as the density in front of the forward shock continu- 
ously increases and the heated matter compensates in the 
flux for the deceleration of the forward shock. This remains 
until the forward shock of the second shell and the BM shell 
merge at time t b s = 0.31 days. At that time the optical 
light curve presents a change in the slope and the flux starts 
decreasing as the emission now originates from the merged 
shell while propagating into the ISM. 

In figure [6] we plot the optical light curve for case 4 
considering a hard-edged jet with opening angle 28h = 2 
during the period of the flare and the emission from several 
snapshots before, after and during the collision of the two 
shells. The emission from each snapshot consists of a BM 
shell contribution which peaks at early observer times and 
the contribution from the second shell which peaks during 
the flare and at later times. A comparison between the light 
curve and the emission snapshots shows that immediately 
after the forward shock of the second shell is created the 
flare is observed (t e = 7.64 x 10 6 s). The stage of the for- 
ward shock deceleration and the flux reduction can be seen 
from t e = 7.64 x 10 6 to t e = 1.01 x 10 7 s. As soon as the 
forward shock separates from the second shell, the shape of 
the emission arising from the forward shock becomes sharper 
(t e = 1.05 x 10 7 s). The third stage of the flare where the 
forward shock propagates inside the BM matter shapes the 
plateau which is observed from t bs = 0.25 to t bs = 0.31 
days and can be seen from t e = 1.05 x 10 7 to t e = 1.26 x 10 7 
s. At that emission time the merger of the two shells has 
almost completed and the flux starts decaying. 



3.4 Time delay observed between optical and 
radio flares. 

Comparing the optical to the radio light curves, we notice 
that the flaring activity occurs with a distinct time delay 
(approximately 0.1 days) for the latter ones (Fig. [7]). The 
reason for this is the ssa mechanism. For optical emission 
the jet is optically thin and the contribution from the sec- 
ond shell is obtained as soon as the forward shock is created 
while propagating in the BM shell. Below the self-absorption 
frequency though the jet behaves differently. In the radio the 
jet is optically thick due to ssa mechanism and the merger 
becomes visible only after the collision has nearly completed. 
This results in observing the plateau at the second stage of 
the optical flare and a sharp peak in the radio flare. This 
highlights the significance of taking into account ssa when 
calculating radio light curves. In fact any variability result- 
ing from changes in the fluid conditions m ay manifest in 
a chromatic fashion (Ivan Eerten et ai1l2010bh . For that rea- 
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Figure 7. Optical and radio light curve comparison for case 4, 
assuming a jet opening angle 28^ = 2. The optical flare precedes 
the radio one by 0.1 days. The three different stages of the flare 
evolution can be seen on the optical flare. 



son the jet break is significantly postponed in the radio light 
curves (Figs. [4] and [5} ■ 



3.5 Emission images 

By directly plotting the relative contributions to the light 
curve from different parts of the fluid, the reasons for the 
differences between radio and optical light curves become 
even more obvious. The flux at a given observer time is 
obtained by solving the linear radiative transfer equation 
through the evolving fluid for a large set of rays. By sepa- 
rately storing the local contributions to the emerging rays, 
we have created emission images showing exactly the rela- 
tive contributions of different parts of the jet. We produce 
emission images for different observer time and opening an- 
gle, for optical (5 x 10 14 Hz) and radio (10 8 Hz) frequencies 
covering the time before, during and after the flaring activ- 
ity is observed. In all the images the jet axis is aligned to 
the horizontal direction and the observer is at the right end 
of the horizontal axis. Which area of the jet is the main con- 
tributing area to the emission is strongly frequency depen- 
dent. When the frequency of observation lies well above the 
self absorption frequency, the system is optically thin and 
the main contribution to the emission is from early emission 
times and from high emission angles on the jet. For lower 
frequencies the system becomes optically thick due to self 
absorption, hence the main contribution to the emission is 
from later emission times and the emitting region shifts to 
lower emission angles closer to the jet axis. This behaviour 
is observed throughout the figures 1 81101 where we plot the 
ring-integrated, absorption-corrected local emission coeffi- 
cients, as a difference in the contrast between optical and 
radio emission. 

In Fig. [S]we plot the optical and radio emission images 
at observer time t b s = 10 days, long after the merger has 
completed, for case 4 for a spherical explosion and for a hard- 
edged jet of opening angle 29^ = 2. We observe that differ- 
ent parts of the merged shell contribute to different parts 
of the spectrum. The optical contribution emerges mainly 
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emission image in case 4 at time = 10 days. The upper pic- 
tures correspond to a spherical explosion and the lower ones to a 
jet with opening angle 28^ = 2. The horizontal axis scales from 
2xl0 17 to 1.5X10 18 cm and the vertical from5xl0 15 to 3.5xl0 16 
cm. The main contribution area for the optical image comes from 
higher emission angles while for the radio image lower emission 
angles contribute the most. For a hard-edged jet (lower graphs) 
the radio image appears stronger, since the main contribution 
area to the optical at the back of the jet is not taken into account. 
The remainder of the second shell significantly heated after the 
traverse of the reverse shock reveals its contribution to the radio 
image. 



from high emission angles of the merged shell whereas ra- 
dio emission comes mainly from lower emission angles close 
to the jet axis. This behaviour is also encountered in the 
emission images derived for the hard-edged jet. In this case, 
high emission angles are excluded from the calculation of 
the emission, which affects the optical more strongly than 
the radio emission. A contribution from the remainder of 
the second shell is observed at the radio emission image of 
the spherical explosion. Although significantly reduced due 
to ssa and almost two orders of magnitude less than the 
emission originating from the merged shell, this contribu- 
tion owns its existence to the traverse of the reverse shock 
(Fig. [8]). As the reverse shock crosses the second shell the 
Lorentz factor decreases, and by time t &s = 10 days, when 
the shell has been significantly decelerated, the peak fre- 
quency of the synchrotron spectrum originating from that 
region is shifted to lower frequencies. The stronger the re- 
verse shock the stronger the deceleration of the shell. The 
contribution of the second shell is therefore expected to be 
higher. 

In Fig. [9] we show the evolution of the two shells before, 
during and after the collision phase for case 2, for both radio 
and optical frequencies assuming opening angles of the jet 
2 and 20 degrees (rather than 5, to better capture high- 
angle emission and illustrate the contrast between narrow 
and wide jets). The collision appears between observer times 
tobs = 0.28 and t bs ~ 0.35 days which verifies that the break 
in the light curves occurs at the same time as the collision 
of the two shells. From the emission images it becomes clear 
that the smearing out of the flare from early emission time 
contribution becomes less the smaller the opening angle of 
the jet is. Comparing the two different cases in figure [9] (case 
2) and figure \W\ (case 4) it is evident that in the optical 
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Figure 9. Optical (5 X 10 14 Hz) and radio (10 8 Hz) emission 
images of the two-shell system during the flaring activity for case 
2 and jet opening angle 2 degrees (top figures) and 20 degrees 
(bottom figures). The horizontal axis scales from 2 X 10 17 to 5 X 
10 17 cm and the vertical from 2 X 10 13 to 8.5 X 10 15 cm for 
the 2 degrees jet and from 2 X 10 13 to 2.5 X 10 16 cm for the 20 
degrees jet. For each subfigure the bottom images correspond to 
tobs = 0.23 days, which is the time the sudden rise of the flux is 
observed, the middle ones to the weak decay at t a i, s = 0.28 days 
corresponding to the propagation of the forward shock into the 
BM medium, and the top ones to the fast decay at t b s = 0.35 
days, once the merger has completed. 
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Figure 10. Optical (5 X 10 14 Hz) and radio (10 8 Hz) emission 
images of the two-shell system during the flaring activity for case 
4 and jet opening angle 2 degrees (top figures) and 20 degrees 
(bottom figures). The horizontal axis scales from 2 X 10 17 to 5 X 
10 17 cm and the vertical from 2 X 10 13 to 8.5 X 10 15 cm for 
the 2 degrees jet and from 2 X 10 13 to 2.5 X 10 16 cm for the 20 
degrees jet. For each subfigure the bottom images correspond to 
t obs = 0.23 days, which is the time the sudden rise of the flux is 
observed, the middle ones to the weak decay at t D i, s = 0.28 days 
corresponding to the propagation of the forward shock into the 
BM medium, and the top ones to the fast decay at t a i, s = 0.35 
days, once the merger has completed. 



case and for large opening angles, early time contribution 
dominates the emission and drowns out the flaring activity. 
For small opening angles this effect is suppressed revealing 
this way the contribution of the second shell while it collides 
with the BM shell. In the radio emission images, where the 
main contributing area to the emission is transferred to lower 
emission angles, we do not observe this behaviour. Instead 
as described in section 3.2, in the radio frequency the jet 
is optically thick due to ssa mechanism and the merger can 
be seen only when the collision has almost completed. This 
difference in behaviour between optical and radio, manifests 
itself as a sharp rise in the radio light curve and a plateau 
at the flare in the optical light curve (see also Fig. [7)1 . 

From the emission images in fig. [9] and for a jet with an 
opening angle 20 degrees we derive the conclusion that for 
higher frequencies (optical, X-rays), for which the main con- 
tributing region to the observed flux shifts to higher emission 
angles on the jet, a double ring shaped image should appear 
in GRB late afterglow observations. In this type of image 
the inner ring carries information from the forward shock 
emission of the second shell, whereas the outer ring is de- 
termined from the emission at earlier stages of the forward 
shell. For radio frequencies where the observed flux arrives 
mainly from emission angles close to the jet axis, the distinct 
contribution of the two shells can only be observed after the 
merger has completed. 



4 DISCUSSION AND CONCLUSIONS 

We performed high resolution numerical simulations of late 
collisions between two ultra relativistic shells and produced 
optical and radio light curves and emission images for a 
spherical explosion case and different opening angles of a 
hard-edged jet. The AMR technique allowed us to reach 
high resolution and properly resolve the shocks developed 
during the merger of the shells. 

The simulations have shown that different values of the 
Lorentz factor and energy of the second shell can signifi- 
cantly change the characteristics of the variability on the 
light curves. We demonstrate that for small opening angles, 
for which the flux is not smeared out by early emission time 
contribution, a flare appears at the light curve. The onset of 
the flare is found to be strongly dependent on the strength 
of the forward and reverse shock of the second shell and to 
become stronger the higher the Lorentz factor and energy of 
the second shell are. For case 4 in which the Lorentz factor 
and energy of the second shell are the highest in our sim- 
ulations, the relative increase of the flux in respect to the 
underlying afterglow is AF/F — 3.95 for the optical and 
AF/F = 1.04 for the radio. The timescale of the variation 
of the flare is AT /T = 1.08 and AT /T = 1.90 respectively. 
The shape of the flare is understood through the dynamical 
simulations and the different stages at the light curve are 
associated with different parts of the collision process. We 
show that the difference in shape between the optical and 
radio flare as well as the time difference observed between 
them is a direct result of the ssa mechanism and the an- 
gle dependence of the emission. We predict that this type 
of behaviour should appear in late afterglow observations 
as a two-ring feature in spatially resolved optical emission 
images, although for the time being the resolution required 
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would be unrealistic. Athough a detailed numerical approach 
is required to fit observational data with the numerical re- 
sults, a straightforward comparison was made between case 
4 of our simulations and the optical flare observed in the 
afterglow of GRB 060206 showing strong resemblance in 
the magnitude and time variation of the flare. New low- 
frequency facilities such as ALMA are expected to provide 
enough data to compare with the radio light curves as well. 
The details which determine the jet collimation as the jet 
propagates into the circumburst medium are currently not 
well understood. In this work we have used the same open- 
ing angle for both shells and we defer analysis of a scenario 
with different opening angles to future work. Such a study 
would require detailed simulations in 2D. 

We have not explicitly discussed X-ray flares, although 
for these the largest amount of observational data is avail- 
able. The X-ray emission is influenced by electron cooling 
(i.e. lies above the cooling break in the synchrotron spec- 
trum) and is therefore also more sensitive to the details of 
particle acceleration than optical and radio emission. Espe- 
cially in the case of multiple shocks, it becomes difficult to 
implement an approach to radiation that correctly captures 
all relevant physics in a simple parametrization. Many after- 
glows show X-ray flares superimposed on a broken power- 
law where a single slope describes the region before and af- 
ter the flare and are often characterized by small timescale 
variations, AT/T < 1 (e.g. Burrows et al. 2007). How- 
ever, there are occasions, such as the GRB 060206, where 
the afterglow light curve presents a clear difference between 
the pre and post flare region with a significantly increased 
timescale variation, indicating an energy injection to the ex- 
ternal sliock_j : s_a_rjr^ablecause for the generation of the 
flare (jMonfardini et al.ll2006r ). 
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